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The process of pattern formation in the two dimensional Swift-Hohenberg equa- 
tion is examined through numerical and analytic methods. Dynamic scaling rela- 
tionships are developed for the collective ordering of convective rolls in the limit of 
infinite aspect ratio. The stationary solutions are shown to be strongly influenced by 
the strength of noise. Stationary states for small and large noise strengths appear to 
be quasi-ordered and disordered respectively. The dynamics of ordering from an ini- 
tially inhomogeneous state is very slow in the former case and fast in the latter. Both 
numerical and analytic calculations indicate that the slow dynamics can be charac- 
terized by a simple scaling relationship, with a characteristic dynamic exponent of 
1/4 in the intermediate time regime. 
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I. INTRODUCTION 



Convective instabilities in liquid systems provide an interesting example of nonlinear, 
nonequilibrium processes. An example of such phenomena that has received considerable 
attention in recent years is the onset and formation of roll patterns in Rayleigh-Benard 



above the convective threshold, for a system of infinite horizontal dimension. Once above the 
convective threshold, a complex pattern of convective rolls emerges. By studying extremely 
large systems, we focus our attention on the collective behavior of a large number of rolls. If 
all sources of fluctuations are small (let F denote their amplitude) the rolls form small locally 
ordered domains (defined by regions of rolls with the same orientation) that later reorient 
to become parallel to rolls of larger domains in their immediate vicinity. Dynamic scaling 
relationships analogous to those found in spinodal decomposition are used to understand 
the transient dynamics of the collective ordering, while ideas based on critical phenomena 
are used to analyze the asymptotic steady states. 

The SH equation was developed by considering the set of equations for a simple fluid 
in the Boussinesq approximation, bounded by two infinite horizontal plates separated by a 
distance rf, at temperatures Ti and Ti + AT respectively. For values of the Rayleigh number 
R larger than a critical Rayleigh number i?c, an instability occurs giving rise to convective 
rolls. The solution of the linearized equations for the velocity field perpendicular to the 
plates Vzi and for the temperature field 9 (or more precisely for the deviations from a linear 
temperature gradient between the plates) contains stable and unstable eigenvalues. The 
SH equation is obtained by neglecting all terms proportional to the stable mode and by 
considering only the wavelengths near the most unstable wavelength obtained in the linear 
analysis. The SH equation is asymptotically correct in the limit R — Rc- In dimensionless 
units, the SH equation reads. 



convection p]-p^. In this paper we consider a simple model of this process, the Swift- 
Hohenberg (SH) equation, to study the situation in which a nonconvective state is brought 




(1) 
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where i/j is a scalar, two-dimensionald field related to the amplitude of the eigenfunction 
corresponding to the unstable mode and is commensurate with the convective rolls. The 
quantity e = {R — Rc) / Rc acts as a control parameter, and r/ is a random field, Gaussianly 
distributed, with zero mean and correlations, 

< r/(f, t)r]{r', t') >= 2F5(f - 7^)6(1 - t'), (2) 

where F is the intensity of the noise. 

The salient feature of the SH equation is that the Lyapunov functional associated with Eq. 
(1) is minimized by a one dimensional periodic function which corresponds to a configuration 
comprised of straight, parallel rolls. The appearance of such "striped" patterns is a common 
feature in nature and has also been observed in magnetic films, block co-polymers, the visual 
cortex, liquid crystals, microemulsions and eutectic growth. Equation (P can be expressed 
in terms of a Lyapunov functional in the following manner, 

where JF is given by, 

jr^y^/f ^ (^e_ (l + V2)^^^/2 + V'V4 . (4) 
As noted in the original work of Swift and Hohenberg this model does not fall within the 



classification scheme of Halperin and Hohenberg |jT5[ . Although other systems select a char- 
acteristic wavelength which is finite (as in order-disorder or anti- ferromagnetic transitions), 
this model differs in that the pattern of rolls is rotationally invariant. Thus small changes 
in orientation are associated with very small changes of JF. Another important feature of 
Eq. is that ip is a non-conserved variable since / df'il){r,t) can vary in time. In Section 
(II) a Lyapunov functional with these two generic properties is considered, in order to both 
study the SH equation and to establish the extent to which our work may apply to other 
systems. 

For the one- dimensional SH equation, Pomeau and Manneville |^ have obtained an 
approximate solution ip^"^ that minimizes T in the limit of small e. This solution will be 



used to estimate the size of various quantities that enter into the analytic calculation given 
in Section (II). The solution is given in a power series in e, and is, 

oo 

^'''(a;) = E«iSin(/eo[2« + l]x), (5) 

i=0 

where, ao(^o) = ^4c<;(A;o)/3, a.i(/co) = — '^o(^o)^/(4ci;(3/co)), and uj{ko) = e — (1 — /cq)^, and 
kf) = 1 — e^/1024. In two dimensions, a solution that minimizes is for example ip{x,y) = 
ip^'^{x). This solution corresponds to a set of straight rolls parallel to the y direction. 
By analogy with a liquid crystal, we will refer to such a configuration as smectic since 
it contains both long-range translational and orientational order. Although this solution 
minimizes JF, it is doubtful that it could be obtained via numerical simulation of the SH 
equation from an initially random state due to the difficulty in removing defects in the 
absence of fluctuations. In fact, Pomeau and Zaleski and Kramer and Zimmermann 
have derived other stationary solutions to the one- dimensional SH equation that are not a 
global minimum of JF and involve a local shift in the phase of ijj (e.g., ijj sin{kox) for x < —a 
and sin{koX + vr) for x > a). 

The presence of fluctuations (i.e., F ^ 0) changes matters significantly. In one dimension 
long range order is broken. A numerical solution of the stochastic SH equation in one 



dimension by Vinals et al. has shown that a finite e' leads to a diffuse peak in S{k,t). 
This one dimensional result has important implications for the two dimensional stationary 
state. If the two dimensional stationary state is a set of parallel rolls, then no long range 
translational order can exist, since the one dimensional solution should be valid normal to 
the rolls. The two dimensional equilibrium state could still, however, exhibit long range 
orientational order. One interesting feature of the SH equation is that the orientational 
order parameter (i.e., a vector that lies normal to the convective rolls) is continuous, similar 
to the spin field in the XY model. Of course the distance between spins in the XY model 
is fixed on a lattice, whereas in the SH equation the distance between convective rolls can 



vary. Since the XY model is the prototypical model of a Kosterlitz-Thouless |]TB[ transition 
it is possible that the SH model contains such a Kosterlitz-Thouless phase. In fact Toner 



and Nelson [|TT] , p!2| have suggested that this could be the case. In Sections (II) and (III) a 
more detailed discussion of the stationary states will be presented. 

The dynamic evolution to the stationary states is realized through the enlargement of 
regions of convective rolls with the same orientation and the elimination of defects. A 
singular perturbation solution ||^ to the SH equation predicts that correlations in ip obey 
the following scaling relationship, 

sik,t) = t^f{[k-ko]n (6) 

where, S{k,t) = J2k < l^(^!^)P > /Z^fc; ^ is the wavevector, a; is a dynamic scaling 
exponent, and the orientation of k have been averaged over. This relationship should be valid 
for extremely large systems, comprised of many domains of different orientation. Equation 
(^) is analogous to the scaling relationships observed in spinodal decomposition and order- 
disorder transitions. In domain growth phenomena the structure factor has been found to 
obey the following dynamic scaling relationship [jl^|T8|, S{k,t) = R{ty f [{kR{t)], where 
R{t) oc is the average domain size and d is the dimension of the system. The differences 
between this result and Eq. (^ are due to the existence of a non-scaling length: the roll 
width. The specific value of x has been shown to be of considerable importance in first order 
phase transitions, and is particular to a given universality class. Systems with a conserved or 
non-conserved scalar order parameter fall into classes characterized by x = 1/3 and x = 1/2 
respectively. 

In Sections (II) and (III) analytic and numerical methods are used to estimate the dy- 
namic exponent x for the SH equation. The analytic calculation is an expansion around a 
curved set of rolls of varying width, and is an extension of the interfacial methods used in 



first order phase transitions. Some of the results of Toner and Nelson [|ri|,|l3 alluded to in 
the preceding paragraph are recovered in Section (II), and are shown to be consistent with 
the numerical results given in Section (III). In the final section a discussion of these results 
is presented and comparisons with other systems is made. 
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II. DYNAMICS OF CURVED ROLLS 



Insight into the ordering dynamics can be obtained by considering the relaxation of a 
set of convective rolls in which the orientation and spacing varies slowly in space. This 
description precludes the explicit consideration of defects, which will be discussed later. 
A typical pattern that falls within this description is shown in Fig. (|l|) which was taken 
directly from the numerical calculations to follow. In the absence of defects, the dynamic 
evolution can roughly be separated into three distinct mechanisms: the relaxation of the 
local curvature of the rolls, the relaxation of the distance between consecutive rolls, and 
the relaxation in the functional form of ip (to be described later). The latter fluctuations 
will be seen to decay much faster than the other mechanisms which describe changes in the 
positions of the roles. To begin the calculation a Lyapunov functional J^{'4'} of a scalar 
function ip of the form, 



= I df[-i; (/ (V^)) i;/2 + giij)] , (7) 



is considered, with the restrictions that JF is minimized (in the absence of noise) by a roll 
structure with characteristic wavenumber ko, is non-conserved and g is an even function 
of ip. Subject to these restrictions the functions / and g can be written, / = Y,i=o o^iiV'^Y 
and g = I]i=i Cj(V'^)*^^- The SH equation can be recovered by choosing = e — l,ai = 
— 2,a2 = —1 and aj>3 = 0, and Ci = 1/4 and Cj>i = 0. We restrict our analysis to even 
powers of ip due to the symmetry of roll patterns, while pointing out that odd powers of ip 
can lead to interesting behavior such as stationary solutions with hexagonal symmetry. The 
dynamical evolution of ip is given by, 

^ = + ^(r, t) = /(V^)^ - Sg/S^ + v{f, t). (8) 

Other systems that are modeled by Eqs. (|^) and (^ include pattern formation in the visual 
cortex [|1^] and a model of uniaxial ferromagnetic films pO| . 

In order to describe the relaxation of the configuration illustrated in Fig. (|T]), the 
Cartesian coordinates x and y are mapped onto curvlinear coordinates u and s as indi- 
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cated in the figure. In this new coordinate system, a vector f{x, y) = xx + yy becomes 
fVn(s, u) = Rm{s) + {um — 2m7i / ko)nm{s) , whcrc Rm{,s) is the location of the m*^ roll, nm{s) 
is the normal to the m*^ roll, Um is restricted to the region (2m — l)7r/A;o < u < (2m + l)7r/fco, 
and m = 0,±1,±2,.... The coordinate system is shown in Fig. (|1|). In the new coordinates 
the Laplacian operator becomes, 

^2 _ _^ ^ ^ ^ 1 {u^ - mTc/ko)Ks d 



dv? {1 + {um — 'rn'K/ko)K) du {1 + {um — rmx /kQ)Ky ds"^ (1 + — mTr/fco)^)'^ 9s 

(9) 

where the curvature k = —d6{s)/ds, Kg = dn/ds and 6{s) is the angle between the tangent 
to the Um = 2m7c/kQ curve and the x axis. Assuming that the curvature of all the individual 
rolls is small, i.e., ttk/Zcq << 1, Eq. @ becomes, 

_o 92 9 92 



9^2 du^ 9^2 ^ 

Coupling between rolls is incorporated by considering the fluctuations in the separation 
between them. The field ip is expanded around the one dimensional stationary solution in 
the following manner, 

ij{r, t) = ^^''{umir, t) + h{um, t)) + 5i,{r, t), (11) 

where, ip^'^ is defined by the equation, 

/(92/9..2)^^'^-^U..=0. (12) 

This decomposition describes the relaxation of the local curvature through the curvature 
dependence of the Laplacian (Eq. ([T0|)), and the relaxation of the separation between rolls 
through the function h. dip takes into account fluctuations about the one dimensional 
functional form of ip described by Eq. (p]2D. The validity of this expansion depends on 
the assumption that the one-dimensional solution ip^'^ is valid normal to the lines defined 
hy ip = 0. The numerical results to follow verify that this is quite a good approximation 
at late times. As can be seen in Fig. (P, tp'^'^ is found to be a good approximation to 



ifj in the direction normal to the rolls for distances of 6 to 20 consecutive wavelengths. 
The assumption that h is small precludes a rapid change in the phase of ip that would be 
associated with defects. 

We next substitute Eqs. (|T0|) and (0) into Eq. (H), and expand to lowest order in k, h, 
6ip, and their derivatives. This gives the result, 

°° °o 4(4 — 1 ^ 

{dij^'^/dw){dw/dt) + {d6i)/dt) = KY,iai{d/dwf'-^i}^'^ + KssY. o, ai{d/dwf'"^i)^'^ 

i=l i=2 ^■ 

oo 

+ dh/duY,{2i)ai{d/dwf'i)^'^ 

i=l 
•t=l 

- \5^T/6ij'') + ^ + • • • (13) 

where w = Um + h{um)- The coefficients of the terms proportional to k and dh/du are iden- 
tically zero if ip^''' is of the form sin(fcou^) or cos(fcow) with defined by {df{k)/dk)\k=ko = 0. 
For the SH equation replacing ip^'^ by a sinusoidal function is an extremely good approx- 
imation for small e as was shown by Pomeau and Manneville 0. The coefficient of the 
term proportional to 6ip is less than zero since tp^"^ is a minimum of JF which implies 

{S^J'/dij^^^^^^ > 0. 

The evolution of 6tp can be obtained by introducing the projection operator, 
Vh = ko/{2n) S^Z-Z/k^ dw. Applying Vh to Eq. (P gives, 

dS^/dt = (/(V^) - VH{S'gW/S^')\^=^..) 5^ + V, (14) 

where, < T]'{n, Si,ti)T]'{m, S2,t2) >= 2F{Sn,m/{'^Tr/ko))S{si - S2)S{ti - 12). In deriving this 
equation, it was explicitly assumed that the length scale over which 5ip varies is much greater 
than the roll wavelength, that is, 5ip is approximately constant over length scales of the order 
27r//co. In addition it was assumed that the coefficient of dh/du was negligible. For the SH 
equation, Eq. ([l^) becomes, 

d5%p/dt = (-2e - (1 + V^)^) 6tp + r/', (15) 
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if Eq. d^) to lowest order in e is used for ip^'^. Thus Sip decays exponentially in time since 
— 2e — (1 — k'^)^ is less than zero for all k. In the next few paragraphs, it will be shown that 
the other mechanisms decay as a power law in time and consequently the relaxation of the 
initial roll pattern considered is controlled by the motion of the rolls. 

These dynamics can be extracted by applying a different projection operator, i.e., 
Vp = kJ{2iT) dwid^'^/dw). Applying Vp to Eq. ® gives, 

dw/dt = biK + b2K,ss + bsd'^h/du^ + ({u, s, t) (16) 

where, 

fei = ^mi^^, (17) 



^3-E«.^^^^^, (19) 



(Ji = fco/(27r) / dwidij^'^/dw)id'ij^'^/dw'), (20) 

J{2m-l)TT/ko 



and < C{m,s,t)C{n,s',t') >= {2F/ai){6,n,n/{27T/ko))6{s - s')6{t - t'). Equation (|T6]) is the 
main result of this calculation. A simple dimensional analysis of Eq. ([T6|) reveals that two 
different mechanisms provide different relaxational rates: The distance between rolls relaxes 
at a rate of t^^/^, while the curvature of the rolls relaxes at a rate of hit^^/'^ + b2t~^^'^. If 
ip^'^ = A sin(A;o'u), then the coefficient bi is identically zero, and the rolls straighten at a rate 
of which is considerably slower than the relaxation of the fluctuations in the distance 
between rolls. Therefore curvature relaxation is the dominant mechanism in the limit bi = 0. 
Corrections to a sinusoidal ip^'^ are generally small close to onset. This leads to a small but 
finite value of bi (e.g., for e = 0.25, bi ^ 1.0 x 10^^). The low value of bi indicates that 
the asymptotic behavior in which the curvature relaxes with a single power law with an 
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exponent of 1/2 should not be expected until very late times (i.e., of the order t = 10^ for 
e = 0.25). 

Equation (|T6|) also describes the local fluctuations in the position of the rolls. If the 
magnitude of the noise intensity (i.e., 2F/ai) is sufficiently strong, i.e., greater that the 
distance between the rolls, the pattern of parallel rolls would be broken. Hence a transition 
to a disordered state will occur when Fkt oc cri/kl (for the SH equation this corresponds 
to Fkt oc e). It should be noted that a is typically proportional to k^. The constant of 
proportionality is difficult to obtain, moreover it would depend on the existence of defects. 



The effect of defects in related models has bee studied by Toner and Nelson |Tl|,|12 
Their analysis begins with a generic free energy of the form, 

G{w) = {B/2) j df {dw/dyf + (dw^/dx^)" . (21) 

Assuming that hi is negligible, Eq. ([I6D can be written in the form, 

dw{x, y, t)/dt ^ -SG{w)/Sw + (. (22) 

where the free energy G{w) is given by Eq. (P), with 5 = 63 and = -fea/^s- Toner 



and Nelson [|n|,0 have shown that defects in conjunction with the generic free energy 
given in Eq. (^) leads to algebraic decay in the orientational correlation function (C(r)) 
, i.e., C(r) =< e2^<^(^^e2*'^(°) >^ \r\-\ where = -dw/dx and ^ v = -2F/{nb2cri). 
The algebraic decay with distance with a F dependent exponent implies that the state is a 



Kosterlitz-Thouless phase. By using the theory of Kosterlitz and Thouless |jT6[ a transition 
to a disordered state occurs when v = 1/4 or Fkt = — vr620"i/8. For the SH equation 
(assuming ip^^ = J Ae/3sm{w)) this relationship becomes, Fkt = vre/12. If an analogy with 



liquid crystals can be made, as suggested by Toner and Nelson |TT|,|T2[, various phases can be 
identified; an isotropic phase with short-range order (for F > Fkt), a nematic or Kosterlitz- 
Thouless phase with "quasi" long-range orientational order (for < F < Fkt) and a smectic 
phase with long-range orientational and translational order (at F = 0). Portions of typical 
final configurations along with the phase diagram are shown in Fig. (|^). An implication 
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of this result is that in order to observe a scahng relationship an orientational correlation 
function should be measured, however S{k,t) should provide an accurate description of 
the ordering dynamics when the translational correlation length is larger than the average 
domain size. 

The results of the preceding paragraphs indicate that the dynamic in the absence of 
defects is controlled by the straightening of rolls, and the rate associated with this mechanism 
is during an intermediate time regime. The effect of defects on the dynamics has been 



considered for a similar model (the XY model) by Kawasaki ||2T| and Loft and DeGrand p2 



Dimensional arguments and numerical simulations of the XY model led Kawasaki ||21| and 
Loft and DeGrand respectively to the conclusion that the separate mechanism of defect 
recombination occurs at a rate of If the annealing away of defects in the SH equation is 

the same as in the XY model then the curvature relaxation mechanism should still dominate 
in the intermediate time regime since it is the slower mechanism. The expected crossover to 
X = 1/2 at late times might however be influenced by defects. The more important influence 
of defects is the breaking of translational order and creation of a disordered state above Fkt- 
The idea of expanding ip around an almost parallel roll configuration has been used by 
others, including Toner and Nelson []T2[ and Ahlers et al. [^. In these works ip{x,y) is ex- 
panded around sin{kQx). Our work differs in that ip is expanded around the one dimensional 
solution of the SH equation (which is not in general a sinusoidal function) in a coordinate 
system commensurate with a curved pattern of rolls. In our work ip^'^ is assumed to hold 
in the normal direction to the line defined hj ip = 0, not simply in the x direction. In 
addition, the projection operator techniques not only decouple fluctuations in the functional 
form of ip from curvature and wavelength relaxation, but also show that the fluctuations in 
the functional form of ip are irrelevant since their decay is extremely fast. 
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III. NUMERICAL SOLUTION 



The numerical results to be presented were obtained by discretizing both space and 
time derivatives in Eq. (|^). Euler's method was used to discretize the time derivative and 
the approximation for the Laplacian included contributions from nearest and next-nearest 
neighbors. Schematically the numerical algorithm can be written, 



(e-[l + VT)^(^,J,n) 



-'ip{i,j,nf + r]{i,j,n) 



(23) 



where, the indices represent the coordinates {x,y) and the index n represents time. 
The Laplacian is evaluated using the following discrete operator, 

V^^(z, j) = ( ^ E +^ E -3) ^P{^,J)/{Ax)^ (24) 

\ (nn) (nnn) / 

where, the notation (nn) and (nnn) refers to the nearest and next nearest neighbors to the 
site respectively. Absolute time and spatial coordinates are recovered by the simple 
relationships, t = nAt and r = {ix + jy)Ax. 

As in most numerical simulations of this nature the choice of (At, Ax) is dictated by the 
conflicting constraints imposed by the need for numerical accuracy and the finite compu- 
tational power available. The former constraint requires (At, Ax) to be vanishingly small, 
while the latter requires the opposite. In practical terms the size of Ax is limited by the 
smallest length scale in the problem. Typically the choice for Ax imposes a restriction on 
At which can be obtained from a linear stability analysis of the discrete map Eq. (pSj). In 



some problems there may also exist non-linear numerical instabilities that further restrict 
the size of At. For the SH equation Ax must be smaller than the wavelength selected, which 
is of the order 2tt. Ax was chosen to be 2^/8 to satisfy this constraint. The linear solution 
of the discrete map given in Eq. (|23D contains a numerical instability (i.e., subharmonic 
bifurcation) that occurs when, 

2(Ax)^ 

^ ^ (4- (Ax)2)2 -e(Ax)4" ^ ^ 
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Thus to avoid numerical instability in the linear analysis for e = 0.25 and Ax = 27r/8, At 
must be less than 0.067 (for the standard nearest neighbor discrete Laplacian operator At 
should be less than 0.014). Although this analysis does not take into account the noise term 
or the non-linear terms, it does provide a minimum requirement. Under these considerations 
a time step of At = 0.05 was chosen and no numerical instabilities were observed. Test runs 
with a smaller Ax indicate that the dynamics of S{k, t) are insensitive to decreasing grid 



size near k = |2^. The bulk of the simulations were run on a 512 x 512 periodic lattice for 
e = 0.25 at F' = 0, 0.05 and 0.075 (where F' = F/ (Ax)^) and averaged over 25 independent 
runs. Smaller systems were used to examine the dynamics (at F' = 0.025,0.065 and 0.1) 
and stationary solutions (at a large number of values F' between and 0.09). For non-zero 
F', ip was initially set to zero, while for F' = 0, ip was chosen to be a random variable 
Gaussianly distributed with zero mean and variance 0.1. 

The results of the numerical simulations suggest that there is a qualitatively different 
behavior between the runs conducted at high and low noise strengths. At high noise strengths 
a steady state was rapidly reached which corresponded to a disordered structure. Figures 
(3a) and (3b) display the configurations at t = 2500 and 10^ respectively for F' = 0.075. The 
striking similarity of these configurations and small domain sizes indicates that a disordered 
stationary state has been achieved. In contrast, at low noise strengths a slow ordering 
of domains was observed for all times probed. Figures (3c), (3d) and (3e) display the 
configurations at t = 10^, 10^ and 10^ respectively for F' = 0.05. A similar time sequence 
is shown for F' = in Figs. (3f), (3g) and (3h). It is important to note that these runs 
had not reached a steady states as the rolls were continuing to order. The pronounced 
difference between the low and high noise strength simulations was also apparent in the 
structure factors and one-point distribution function p{ip). In Fig. (4a) S{k,t = 700) is 
shown for F' = 0.075. This structure factor was statistically indistinguishable from any 
S{k,t) past t = 200. The solid line in this figure is a fit to A/{B + (A;^ — k'^Y), which is 
suggestive of a disordered phase. The structure factors shown in Fig. (4b) and (4c) for 
F' = 0.05 and respectively are significantly different for two reason. First their shape is 
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much sharper, and secondly they are continuing to evolve in time. Further indication of the 
qualitative differences is shown in Fig. (5) which compares p{ip) at the latest times for the 
three noise strengths. The low noise strength one-point distributions are bimodal, while the 
high noise strength p is peaked ai ip = 0. To obtain an estimate of where the transition 
occurs at for e = 0.25, a stationary solution at F' = 0.09 was obtained and then F' was 
decremented in steps of 0.005 at time intervals of 10^. This test reveals a transition in the 
region F^j, = 0.065 — 0.070 which is signaled by a sharp decrease in the peak height of S{k), 
and a crossover from a bimodal to unimodal distribution in p. 

Although the numerical picture is consistent with the predictions given in Section (II), 
it is important to note the limitations of the numerical results. The predictions of smectic, 
nematic and disordered states were not unambiguously verified. What is apparent is that 
the runs above F'j^j, are disordered states. The numerical results also strongly suggest that 
the states below F'j^j^ are qualitatively different from those above. Indeed, the structure 
factors for the smectic and nematic simulations were indistinguishable at the latest times. 
If the F' = 0.05 is a nematic state, one should observe only quasi-long range orientation 
order. In contrast the numerical results indicated growth of translational order up until 
the latest times probed. Thus the numerical results indicate a transition from a disordered 
state to one with long or quasi-long range translational order. Presumably, if later times 
were probed an orientational order parameter would have to be introduced to observe the 
subsequent ordering of domains. At the very least the transition observed corresponds to 
one in which there is a dramatic change in translational order at F^rp ^ 0.0675. 

The dynamics of the collective ordering of rolls was analyzed in terms of the dynamic 
scaling relationship given in Eq. (|]) for F' < F^j,. In Figs. (^) and (Qb) S{k, t) is 
displayed for F' = and 0.05 respectively at several times. The ordinate in these plots is 
S{k,t)/t^ and the abscissa is (/cq — k)t^. If the scaling predictions are correct, the scaled 
form of S{k,t) should overlap for all times, near k = k^. For F' = (see Fig. (|^a)) the 
structure factor scales remarkable well over the time range t = 10^ to 10^, with a growth 
exponent of x = 1/5. Figure (§b) shows that the dynamic scaling hypothesis also works 
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very well for F' = 0.05, if a growth exponent of 1/4 is used. If the scaling hypothesis is 
valid, all length scales (except the convective roll width) should scale with the same dynamic 
exponent. To investigate this hypothesis the height {A), width {w) and various moments 
mn{t) of S{k,t) were calculated. The height and width were determined by fitting the top 
portion of S{k,t) to a Gaussian of the form Ae^^^'^^^'^o)/^)'^) and the moments were defined 
to be mn{t) = fl^^^^^l^f) dk\k — ko\'^S{k,t). In fitting S{k,t), ko was found to be very close 
to 1 which is the value selected by linear theory. The scaling relationship given in Eq. 
implies that oc t""^, w oc and A oc t^. The results are shown in Figs. (^) and 

(^) for F' = and F' = 0.05 respectively. For F' = the exponents measured from all 
lengths scales are within two percent of 1/5, while the exponents are all within one percent 
of 1/4 for F' = 0.05. Although the statistics collected for e = 0.025 and 0.065 were not 
sufficient to obtain accurate values of x, the measured exponents were significantly closer to 
1/4 than 1/5. 

IV. DISCUSSION AND SUMMARY 

The results of the numerical and analytic work provide several predictions for the col- 
lective ordering of rolls in the limit of infinite aspect ratio. The stationary or steady states 
should be strongly infiuenced by any random noise source, such that above some critical 
noise strength (i.e., Fkt) an isotropic state emerges, while below Fkt an nematic or smectic 
state arises. It should be emphasised that the numerical work indicates that the transition is 
characterized by a drastic increase in translational order, but was not sufficient to determine 
the precise nature of the states below Fkt- An independent estimate of the transition line 
{Fkt oc e) was obtained in Section (II). Further investigation at other values of e would be 
useful to establish this relationship. 

The dynamical behavior of the collective ordering below Fkt was found to obey the 
dynamic scaling relationship given in Eq. (|^) for at least three decades in time. In addition, 
a growth exponent of x = 1/4 was observed in the nematic region. This value of x disagrees 
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with an earlier theoretical calculation and we have presented a possible resolution of the 
discrepancy. We are, however, unable to explain the smaller exponent of x = 1/5 observed 
for e = 0. Although this work focused on the stochastic Swift-Hohenberg equation, we expect 
that both features discussed, namely the existence of a transition between a quasi-ordered 
to a disordered state, and the asymptotic dynamical scaling behavior, are generic features 
of two-dimensional systems defined by Eq. (H). 

The difficulty in observing these phenomena in Rayleigh-Benard experiments is that 
previous experiments were conducted on very small systems in which boundary effects play 
an important role. In addition, estimates of the noise strength for real experiments is 
typically much to small to observe the transition. However recent experimental studies of 
electrohydrodynamic convection in nematic liquid crystals [|2^ have been able to detect and 
quantify the amplitude of the fluctuations before onset (fluctuations in simple liquids are 
believed to be too small to produce observable effects such as the ones presented here). 
Furthermore, since typical roll widths in these systems is of the order of microns, large 
aspect ratio samples are readily available. Experimental verification of the issues discussed 
in this paper seems therefore possible. 
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FIGURES 



FIG. 1. 

Curvlinear coordinate system. The points on this figure are zeros oi t/j a,t t — 10^ for 
e = 0.25 and F' — 0.05. The thick sohd fine was fit to the a fine defined by = 0. The thin 
fines were calculated by the equation f = R + rmr/kofi, were, R is the position of the thick 
solid line, fi is the normal to the solid line and m = ±1, ±2, • • •. 

FIG. 2. 

Portions of size 100 x 100 of typical configurations obtained. The large and small 
circles correspond to the large and small (see text) numerical simulations conducted at 
F' = 0.075,0.05 and 0, and F' = 0.025,0.065 and 0.1 respectively. The configurations la- 
beled isotropic, nematic and smectic correspond to F' — 0.075, 0.05 and respectively. In 
all these plots the fines drawn are the lines oi ip{r) — 0. 

FIG. 3. 

Sample configuration are displayed for e = 0.25 at various times and noise strengths. 
In all figures the points correspond to ^ > 0. The configurations shown in (a) and (b) 
correspond to t — 2500 and 10"^ respectively for F' — 0.075. In (c), (d) and (e), the 
configurations correspond to t — 10^, 10^ and 10^ respectively for F' — 0.05. A similar time 
sequence is shown in Figs, (f), (g) and (h) for F' = 0. 

FIG. 4. 

In Fig. (a), S{k,t) is displayed for F' = 0.075 at t = 700. In this figure the points 
correspond to the numerical result and the sofid fine is a fit to A/ {B + (/c^ — /cq)^). In Figs, 
(b) and (c) S{k,t) is shown at F' — 0.05 and 0.0 respectively. In both these figures, the 
curves (from bottom to top at /c = 1) correspond to t = 10^, 10'^ and lO''. 
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FIG. 5. 



A comparison of the one point distribution function for F' — 0, 0.05 and 0.075 at e = 0.25. 
The sohd, dotted and dashed hues correspond to F' = at t = 10^, F' — 0.05 t — 10^ 
and F' = 0.075 at t = 700 respectively. 

FIG. 6. 

The dynamic scahng of S{k, t) is shown for F' = and 0.05 in Figs, (a) and (b) 
respectively. In both figures the open squares, crosses and solid circles correspond to t — 
10^ 10^ and 10^ respectively. 

FIG. 7. 

The amplitude A{t), width w{t) and first five moments mn{t) of S{k, t) as defined in the 
text are displayed for F' = and F' = 0.05 in Figs, (a) and (b) respectively. The solid lines 
are presented as guides to the eye and all have slope x — 1/5 in Fig. (a) and x = 1/4 in 
Fig. (b). 
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